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We  investigated  multilayer  and  bilayer  Ni/Si  thin  hlms  by  nanodifferential  scanning  calorimetry 
(nano-DSC)  at  ultrarapid  scan  rates,  in  a  temperature-time  regime  not  accessible  with  conventional 
apparatus.  DSC  experiments  were  completed  at  slower  scan  rates  as  well,  where  it  was  possible  to 
conduct  parallel  rapid  thermal  annealing  experiments  for  comparison.  Postexperimental 
characterization  was  accomplished  by  x-ray  diffraction,  and  by  transmission  electron  microscopy 
(TEM)  and  energy-hltered  TEM  of  thin  cross  sections  prepared  by  focused  ion  beam  milling.  We 
found  that  rate  of  heating  has  a  profound  effect  on  the  resulting  micro  structure,  as  well  as  on  the 
DSC  signal.  After  heating  to  560  °C  at  120  °C/s,  the  general  microstructure  of  the  multilayer  was 
preserved,  in  spite  of  extensive  interdiffusion  of  Ni  and  Si.  By  contrast,  after  heating  to  560  °C  at 
16  000  °C/s,  the  multilayer  hlms  were  completely  homogeneous  with  no  evidence  of  the  original 
multilayer  microstructure.  Eor  the  slower  scan  rates,  we  interpret  the  results  as  indicating  a  solid 
state  diffusion-nucleation-growth  process.  At  the  higher  scan  rates,  we  suggest  that  the  temperature 
increased  so  rapidly  that  a  metastable  liquid  was  hrst  formed,  resulting  in  complete  intermixing  of 
the  multilayer,  followed  by  crystallization  to  form  solid  phases.  The  integrated  DSC  enthalpies  for 
both  multilayer  and  bilayer  hlms  are  consistent  with  this  interpretation,  which  is  further  supported 
by  thermodynamic  predictions  of  metastable  Ni/Si  melting  and  solid  state  Ni/Si  interdiffusion.  Our 
results  suggest  that  use  of  heating  rates  >10  000  °C/s  may  open  new  avenues  for  intermetahic 
micro-  and  nanofabrication,  at  temperatures  well  below  those  prevailing  during  explosive 
silicidation.  ©  2009  American  Institute  of  Physics,  [doi:  10. 1063/1. 3254225] 


I.  INTRODUCTION 

As  the  technology  for  complementary  metal  oxide  semi¬ 
conductor  (CMOS)  devices  advances  to  feature  sizes  of  28 
nm  and  below,  replacement  of  the  poly  silicon /Si02  gate 
stack  combination  represents  a  key  paradigm  shift.  The  in¬ 
dustry  is  moving  toward  a  metal  gate/high- ac  dielectric  com¬ 
bination  in  which  fully  silicided  metal  gates  are  receiving  a 
lot  of  attention  as  candidates  for  the  metal  layer.  An  im¬ 
portant  example  is  the  Ni-Si  system  which  consists  of  a 
number  of  phases  with  optimal  electrical  conductivity  and 
work  function.^  In  developing  this  technology,  the  thermal 
stability,  and  thermal  compatibility  with  other  CMOS  pro¬ 
cesses  are  important  factors.^  The  use  of  ultrarapid  thermal 
annealing  processes  which  occur  on  the  millisecond  time 
scale  are  being  investigated  for  the  control  of  dopant  distri¬ 
butions  in  ultrashahow  junctions. It  has  also  recently  been 
demonstrated  that  the  Ni-Si  reaction  can  occur  using  laser 
annealing^’ and  such  ultrafast  processing  of  Ni-Si  repre¬ 
sents  an  important  avenue  for  continuing  Moore’s  law 
progression.^^  While  it  seems  clear  that  the  processes  of  dif¬ 
fusion  and  reaction  that  are  responsible  for  forming  the  sili- 
cides  will  have  different  kinetics,  the  effect  of  ramp  rate  on 
the  reactions  remains  unexplored. 

The  Ni-Si  reaction  has  been  previously  studied  with  Ni 
deposited  on  single  crystal  Si,^  polysilicon,^  and  for  cosput¬ 
tered  Ni  and  Si."^  During  reaction,  more  than  one  phase  may 
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be  present,  for  example,  in  the  work  of  Rivero  et  al.,^  the 
reaction  proceeded  via  a  mixture  of  Ni2Si  and  Ni3Si2  for  a  13 
nm  Ni  patch  deposited  on  a  single  crystal.  In  experiments 
with  Ni  on  polysilicon,^  it  was  found  that  the  ratio  of  the 
thicknesses  of  deposited  Ni  and  Si  determined  the  composi¬ 
tion  of  the  silicide,  with  intermediate  silicides  of  NiSi  and 
Ni3iSii2  formed  during  the  anneal.  Cosputtering  the  two  ma¬ 
terials  also  allowed  control  of  composition  and  was  shown  to 
affect  the  work  function."^  Calorimetric  measurements  have 
proven  useful  at  elucidating  the  processes  during  a  Ni-Si 
reaction,^^  as  have  detailed  kinetic  studies. 

It  was  the  aim  of  this  study  to  explore  the  Ni-Si  reaction 
at  rapid  heating  rates  using  a  combinatorial  approach  with  a 
microelectromechanical  device  that  allowed  for  simultaneous 
heating  and  differential  scanning  calorimetry  (DSC)  mea¬ 
surements.  By  comparing  postreaction  characterization  re¬ 
sults  for  a  bilayer  with  a  multilayer  of  Ni-Si,  each  with  the 
same  total  thickness,  we  were  able  to  investigate  the  role  of 
interdiffusion.  As  expected,  the  reaction  was  enhanced  in  the 
multilayer  sample  due  to  a  short  diffusion  distance.  The  sur¬ 
prising  hnding  is  that  there  is  a  strong  enhancement  of  Ni/Si 
mixing  with  rapid  ramp  rate. 

II.  EXPERIMENTAL  PROCEDURE 

The  nanocalorimeter  shown  in  Eig.  1  acts  like  a  minia¬ 
ture  differential  scanning  calorimeter.^^  The  device  consists 
of  a  66  X  240  X  3.3  /uLuf  thick  inner  suspended  platform  di¬ 
vided  into  sample  and  reference  zones.  The  device  was  fab- 
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FIG.  1.  (Color  online)  Photomicrograph  and  schematic  of  nanocalorimeter. 


ricated  in  an  array  of  six  elements.  Films  are  deposited  on 
the  sample  side  using  a  shadow  mask.  Each  zone  has  an 
integral  poly- Si  heater  and  an  Al/poly-Si  thermopile  for  mea¬ 
suring  the  temperature  difference  between  the  two  zones.  To 
enhance  the  thermal  isolation  of  the  inner  platform,  the  leads 
to  the  thermopile  and  heater  run  along  a  suspended  platform 
that  surrounds  the  inner  platform.  Chips  are  mounted  in 
packages  and  placed  in  a  high  vacuum  system  for  calorim¬ 
etry  scans.  The  thermopile  produces  a  signal  of  1.8  mV/K 
that  yields  1  mK  resolution  at  a  10  kHz  bandwidth.  During  a 
scan,  separate  heater  voltage  (V^)  ramps  were  applied  to  the 
sample  and  reference  heaters,  and  the  thermopile  signal  was 
recorded  using  a  preamplifier  and  digital  oscilloscope.  A  set 
of  devices  coated  with  melting  point  standards  (In,  Sn,  Pb, 
Al-Ge  eutectic)  was  produced  to  calibrate  temperature  and 
enthalpy  for  each  of  the  ramp  rates  used  in  the  Ni-Si  mea¬ 
surements.  DSC  scans  consisted  of  a  ramp  up  to  a  maximum 
temperature,  followed  immediately  by  a  ramp  down  to  room 
temperature.  This  sequence  was  repeated  until  there  was  no 
change  in  the  signal  from  one  scan  to  the  next.  The  signal 
from  the  last  scan  (which  was  due  to  differences  in  the  heat 
capacity  of  the  reacted  sample  and  the  reference  side)  was 
used  as  a  baseline,  which  was  subtracted  from  the  earlier 
scans,  to  provide  difference  scans  showing  the  progression  of 
the  Ni/Si  reaction. 

Films  of  Si  and  Ni  were  deposited  on  the  nanocalorim¬ 
eters  by  magnetron  sputtering  at  a  base  pressure  of  2.7 
X  10“^  Pa.  Two  sets  of  samples  were  prepared  with  the  total 
Ni  and  Si  content  the  same  for  the  two  sets.  Bilayer  samples 
consisted  of  layers  of  100  nm  Ni  and  100  nm  Si,  and  the 
multilayer  films  were  composed  of  ten  sequences  of  (10  nm 
Ni/10  nm  Si).  Films  on  test  coupons  of  Si  wafers  coated  with 
500  nm  of  thermal  oxide  were  prepared  simultaneously  for 
use  in  parallel  RTA  experiments.  Postexperimental  character¬ 
ization  utilized  x-ray  diffraction  (XRD),  focused  ion-beam 
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FIG.  2.  DSC  signal  and  heater  voltage  vs  time  for  2000  °C/s  ramp  rate  for 
the  multilayer  sample. 


milling  and  transmission  electron  microscopy  (FIB/TEM), 
scanning  electron  microscopy  with  energy  dispersive  x-ray 
analysis,  and  reflected  light  microscopy. 

Samples  for  TEM  from  individual  devices  were  cut  and 
thinned  using  FIB  methods.  Compositional  mapping  was  ac¬ 
complished  using  energy-filtered  TEM  (EFTEM)  in  a  JEOL 
3010  ultra-high  resolution  (UHR)  TEM  (300  kV)  equipped 
with  a  Gatan  imaging  filter.  The  Si  L2  3  edge,  Ni  M2  3,  and 
Ni  L23  absorption  edges  were  used.  A  three-window  method 
was  used  to  isolate  the  signal  from  the  background.  Partial 
overlap  of  the  Ni  M2, 3  (65  eV)  and  Si  (100  eV)  edges  pre¬ 
cludes  accurate  background  subtraction  under  the  Si  L2,3 
edge;  however,  semiquantitative  analyses  are  still  possible. 

III.  RESULTS 

A.  Differential  scanning  caiorimetry 

A  representative  DSC  scan  of  a  multilayer  sample  is 
shown  in  Fig.  2,  where  we  plot  a  difference  scan  for  an 
intermediate  heating  rate  of  2000  °C/s.  The  main  features  of 
this  scan  are  a  relatively  large  exothermic  peak  (referred  to 
as  peak  1)  initiating  at  a  heater  voltage  (V^)  of  2.55  V  (cor¬ 
responding  to  234  °C),  followed  by  a  smaller  exothermic 
effect  (peak  2)  starting  at  4.28  V  (483  °C).  We  interpret 
these  effects  as  originating  from  Ni/Si  reactions  activated  at 
different  temperatures.  The  later  exothermic  effect  continues 
beyond  the  maximum  heater  voltage  into  the  cooling  part  of 
the  cycle.  These  effects  are  greatly  diminished  in  scans  be¬ 
yond  the  first,  indicating  that  reaction  was  nearly  complete 
after  the  first  scan. 

In  Fig.  3,  we  show  a  series  of  five  DSC  scans  at  different 
rates.  A  normalization  was  used  for  each  scan  so  that  the 
signals  may  be  conveniently  compared,  as  the  raw  signal 
versus  time  traces  differ  significantly  both  in  duration  and 
amplitude.  The  signal  amplitudes  for  each  scan  rate  were 
normalized  by  the  respective  scan  rate,  and  the  abscissa 
shows  time  divided  by  total  scan  time.  This  ratio  is  the  frac¬ 
tion  of  scan  completed.  Also  shown  is  where 

Vhmax  is  the  maximum  heater  voltage.  The  horizontal  offset 
of  each  DSC  scan  is  adjusted  so  that  all  Vh/V^max  plots 
coincide.  Figure  3  shows  an  interesting  evolution  in  the  fine 
structure  of  the  DSC  peaks  as  the  scan  rate  is  increased.  At 
490  °C/s,  peak  1  shows  double  maxima,  the  higher  voltage 
(higher  temperature)  one  of  which  gradually  diminishes  with 
increasing  scan  rate,  until  no  longer  distinguishable  at 
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FIG.  3.  (Color  online)  DSC  signal  divided  by  the  scan  rate  for  five  scan 
rates  plotted  against  normalized  scan  time.  Also  shown  is  the  normalized 
applied  heater  voltage  V^. 


7950  °C/s.  As  discussed  below,  this  feature  may  indicate 
added  complexity  in  the  reaction  represented  by  peak  1. 

By  integrating  the  areas  under  the  two  peaks  defined  in 
Fig.  2,  it  is  possible  to  get  an  estimate  of  the  enthalpy  of 
reaction  as  a  function  of  scan  rate  by  reference  to  the  Al/Ge 
eutectic  data  mentioned  above.  Results  are  shown  in  Fig.  4, 
where  the  inset  indicates  the  variation  of  the  areas  for  peak  1 
and  peak  2  as  a  function  of  scan  rate.  Note  that  peak  1 
decreases  from  its  maximum  at  the  slowest  scan  rate  to 
nearly  zero  at  the  highest  scan  rate,  whereas  peak  2  shows 
the  reverse  trend.  The  sums  of  the  two  peak  areas  show  a 
monotonic  increase  with  scan  rate,  suggesting  either  that  en¬ 
thalpy  production  and  completion  of  reaction  are  enhanced  at 
high  scan  rates,  or  that  differential  heat  loss  is  greater  at  the 
slow  scan  rates.  For  comparison  we  also  show  in  Fig.  4  the 
cumulative  results  of  similar  experiments  on  a  Ni/Si  bilayer. 
Although  the  data  are  less  extensive,  a  similar  trend  is  evi¬ 
dent  for  the  bilayer. 


FIG.  5.  TEM  images  from  Ni/Si  multilayer  sample  annealed  with  a  ramp 
rate  of  (a)  16  000  °C/s  and  [(b)-(e)]  120  °C/s.  (c)  Magnified  brightfield 
image,  [(d)  and  (e)]  EFTEM  maps  for  Si  and  Ni  from  the  same  area. 


B.  Postexperimental  characterization 

In  Figs.  5(a)  and  5(b)  we  show  microstructures  of  the 
rapidly  scanned  versus  slowly  scanned  multilayered  films,  as 
obtained  by  TEM  bright  field  imaging.  A  large  difference  is 
obvious,  in  that  the  rapidly  scanned  material  [Fig.  5(a)]  is 
essentially  homogeneous,  whereas  the  slow  scan  [Fig.  5(b)] 
retains  a  multilayered  structure.  This  is  further  borne  out  by 
the  energy-filtered  (EFTEM)  compositional  maps  of  Eig. 
5(d)  (Si)  and  Eig.  5(e)  (Ni),  which  show  a  layered  hetero¬ 
structure  of  alternating  Si-  and  Ni-rich  layers.  Detailed 
analysis  of  the  EETEM  spectra  for  the  slow  scan  showed  that 
substantial  interdiffusion  of  Ni  and  Si  occurred,  resulting  in 
relatively  homogeneous  bilayers  of  stoichiometrically  dis¬ 
tinct  Ni-silicides.  In  Eigs.  5(b)  and  5(c),  thin,  discontinuous 
areas  with  high  signal  (bright)  are  evident  at  every  other 
Ni/Si  junction.  These  small  areas  have  the  highest  Si  concen¬ 
tration  and  are  comprised  of  a  third.  Si-rich  phase.  By  con¬ 
trast  with  the  slowly  scanned  material  of  Eig.  5(b),  EETEM 
data  for  the  rapidly  scanned  sample  in  Eig.  5(a)  were  feature¬ 
less  and  showed  no  discernible  compositional  layering;  the 
sample  appeared  to  have  become  homogeneous  with  regard 
to  Ni  and  Si  distribution.  Due  to  the  partial  overlap  of  Ni  and 
Si  edges,  quantitative  information  on  phase  compositions 
was  not  obtained  for  the  slowly  scanned  material. 

Information  on  identity  of  crystalline  phases  in  the  re¬ 
acted  Ni/Si  multilayers  was  obtained  from  two  sources:  (1) 
selected  area  electron  diffraction  (SAED)  data  taken  from 
EIB  samples  of  the  DSC  sensors  and  (2)  XRD  data  on  mul¬ 
tilayered  coupons  subjected  to  a  rapid  thermal  anneal  at 
50  °C/s.  Data  from  the  two  sources  are  consistent  in  that  the 
principal  phases  identified  are  Ni3Si  (major),  (Ni,Si)  solid 
solution,  NiSi2  (major),  and  Si.  No  measurable  difference 
was  observed  between  the  SAED  patterns  from  sections  of 
the  slow  scanned  and  the  rapidly  scanned  DSC  samples. 

IV.  DISCUSSION 


FIG.  4.  (Color  online)  Sum  of  peak  area  (enthalpy)  vs  scan  rate  for 
multilayer  (circles)  and  bilayer  (solid  triangles).  Inset:  peak  1  enthalpy  (dia¬ 
monds),  peak  2  enthalpy  (squares),  and  sum  of  peak  areas  (circles)  vs  scan 
rate  for  multilayer  sample. 


A.  Generalized  interpretation 

Based  on  examination  of  the  combined  TEM,  EETEM, 
SAED,  XRD,  and  DSC  data,  we  developed  working  hypoth- 
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eses  for  events  occurring  during  the  DSC  ramps  as  follows. 
As  the  Ni/Si  multilayers  were  heated,  a  temperature 
(—235  °C,  see  Fig.  2)  was  reached  where  interfacial  reaction 
was  initiated,  resulting,  as  the  temperature  increased  further, 
in  the  first  exothermic  energy  release,  followed  by  a  second 
energy  release  peaking  near  the  maximum  ramp  temperature 
of  560  °C.  As  detailed  in  Fig.  4,  the  intensity  of  these  peaks 
varied  inversely  with  increasing  scan  rate,  such  that  the  ini¬ 
tial  peak  (peak  1),  which  was  not  detectable  at  the  lowest 
scan  rate,  became  the  most  intense  peak  at  the  highest  scan 
rate.  The  second  peak  (peak  2)  which  was  the  most  intense 
peak  at  the  slowest  scan  rate,  became  undetectable  at  the 
highest  scan  rate.  We  associate  peak  1  with  exothermic  Ni/Si 
mixing  (rapid  diffusion  and/or  melting)  at  the  thermody¬ 
namically  unstable  Ni/Si  interface,  and  peak  2  with  second¬ 
ary  crystallization  of  Ni  silicide  phases.  The  competition  be¬ 
tween  rapid  diffusion/melting  and  crystallization/growth  is 
influenced  by  the  scan  rate,  such  that  the  rapid  process  is 
more  pronounced  in  the  rapidly  scanned  materials,  due  to 
kinetic  limitations  on  the  crystallization  rate  of  the  secondary 
solid  phases.  If  the  rapid  mixing  process  results  from  melt¬ 
ing,  it  may  be  to  some  extent  self-sustaining^^  at  the  highest 
scan  rates  where  there  is  less  thermal  loss,  since  for  the  un¬ 
stable  Ni/Si  interface,  the  exothermic  energy  associated  with 
the  negative  enthalpy  of  mixing  of  Ni  and  Si  greatly  exceeds 
the  endothermic  effect  associated  with  the  disorder  of  melt¬ 
ing.  In  the  most  rapidly  scanned  materials  [Fig.  5(a)],  com¬ 
plete  intermixing  of  Ni  and  Si  occurred,  apparently  with 
nearly  simultaneous  crystallization,  since  only  one  discern¬ 
ible  exothermic  peak  was  observed.  The  fundamental  ques¬ 
tion  relates  to  the  origin  of  the  rapid  mixing  process:  melting 
versus  rapid  diffusion  by  a  short  circuit  mechanism  such  as 
grain  boundary  diffusion. 

B.  Metastable  melting 

The  possibility  of  metastable  melting  in  the  rapidly 
scanned  sample  can  be  further  evaluated  by  means  of  a  ther¬ 
modynamic  analysis  of  the  experimental  Ni-Si  phase  dia¬ 
gram.  For  this  purpose,  we  used  THERMOCALC for  the  the- 
modynamic  calculations.  The  calculations  were  completed 
using  the  thermodynamic  descriptions  of  Tokunaga  et  al.}^ 
in  which  the  Gibbs  energies  of  the  disordered  phases  (liquid, 
fee  Ni  solid  solution)  were  modeled  as  regular  solutions  with 
Redlich-Kister  interaction  parameters  fitted  to  the  experi¬ 
mental  data.  For  the  ordered  fee  Ni  phase,  an  additional  con¬ 
tribution  to  the  Gibbs  energy  associated  with  ordering  was 
modeled  with  two  sublattices.  For  the  stoichiometric  phases 
appearing  in  the  Ni-Si  phase  diagram  (including  Ni3Si  and 
NiSi2),  Gibbs  energies  were  calculated  from  the  weighted 
energies  of  the  pure  end  members  combined  with  an  empiri¬ 
cal  term  to  fit  the  experimental  data.  The  resulting  equations 
of  state  allowed  the  calculation  of  metastable  equilibria 
among  any  combination  of  Ni-Si  phases  over  a  wide  range 
of  temperatures. 

In  Fig.  6(a),  the  complete  Ni-Si  diagram  calculated  with 
the  optimized  parameters  of  Tokunaga  et  al  is  presented; 
the  diagram  contains  11  distinct  phases.  In  Fig.  6(b)  we  show 
the  results  of  calculating  metastable  equilibria  among  just 


three  of  the  phases:  the  pure  Ni  and  Si  end  members,  and  the 
liquid  phase.  This  scenario  corresponds  to  the  situation  of 
ultrarapid  heating  where  there  is  no  time  for  growth  of  inter¬ 
mediate  crystalline  phases,  only  formation  of  a  melt.  The 
phase  diagram  of  Fig.  6(b)  indicates  a  deep  trough  of  meta¬ 
stable  liquid  extending  to  low  temperatures,  and  centered  at 
approximate  composition  0.33  mol  fraction  Si.  This  meta¬ 
stable  liquid  encompasses  (fortuitously)  the  bulk  composi¬ 
tion  chosen  for  the  multilayer  films  (0.35  mol  fraction  Si), 
giving  credence  to  the  possibility  of  extensive  metastable 
melting  for  the  rapidly  heated  multilayers.  In  Fig.  6(c)  we 
calculated  metastable  equilibria  with  one  additional  phase: 
the  ordered  fee  solid  solution.  The  presence  of  the  ordered 
solid  solution  phase  significantly  reduces  the  region  available 
for  metastable  melting  and  the  minimum  melting  tempera¬ 
ture  is  raised  to  400  °C. 

C.  Solid  state  diffusion 

If  only  solid  state  diffusion  is  considered  without  melt¬ 
ing,  the  complete  homogenization  of  the  rapidly  heated  mul¬ 
tilayers  relative  to  the  remnant  layering  of  the  slowly  heated 
samples  is  counterintuitive.  Nonetheless,  it  is  useful  to  con¬ 
sider  diffusion  from  a  more  quantitative  perspective.  Assum¬ 
ing  a  composition-independent  diffusivity,  the  Fickian  diffu¬ 
sion  formalism 

dC  d^C 

- =D - y, 

dt  dx^ 

where  C  is  the  concentration  of  Ni,  v  is  the  distance,  and  t  is 
time,  with  boundary  conditions 

dC 

v  =  0  and  x  =  L\ —  =  0, 
dx 

and  initial  conditions 
t  =  0:C  =  kCo, 
where 

k  =  0(x<L/2),  k=\{x>m) 

can  be  solved  numerically  using  an  error  function  solution  to 
give  concentration  profiles  approximating  the  expected  Ni 
and  Si  distribution  at  560  °C  as  a  function  of  time.  For  these 
calculations  we  assumed  that  the  Ni  diffusion  relative  to  Si  is 
the  more  rapid  process.  Since  the  sputtered  films  are  granu¬ 
lar,  we  also  assumed  that  grain  boundary  diffusion  is  the 
dominant  mass  transport  mechanism.  During  reactive  diffu¬ 
sion,  the  diffusivities  of  Ni  through  the  product  phases  must 
be  considered.  Literature  data  are  few,  but  representative 
grain  boundary  diffusivity  [D(gb)]  values  in  different  matri¬ 
ces  are  summarized  in  Table  I,  where  it  can  be  seen  that 
values  vary  widely. 

Specifically,  the  difference  between  Ni  transport  in  NiSi2 
versus  Ni2Si  is  about  six  orders  of  magnitude,  which  would 
have  an  extreme  effect  on  the  calculated  concentration  pro¬ 
files  at  the  Ni/Si  interface.  In  Fig.  7,  we  show  a  diffusion 
profile  calculated  with  D(gb)  for  Ni/NiSi2.  The  Ni  mass 
transport  is  represented  in  the  figure  by  C/Cq,  the  ratio  of  the 
Ni  concentration  after  a  scan  to  the  concentration  before  the 
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Mole  Fraction  Si 


FIG.  6.  (Color  online)  (a)  Calculated  Ni-Si  equilibrium  phase  diagram,  (b)  Diagram  assuming  metastable  equilibria  among  the  pure  Ni  and  Si  end  members 
and  the  liquid  phase,  (c)  Metastable  equilibria  among  Ni  solid  solution,  Si,  ordered  fee  phase,  and  liquid,  (d)  Metastable  equilibria  among  pure  Ni  and  Si  end 
members,  Ni3Si,  NiSi2,  and  liquid. 


scan  shown  as  a  function  of  the  distance  from  the  interface. 
The  calculation  shown  by  the  solid  line  in  Fig.  7  corresponds 
to  the  approximate  time  at  560  °C  which  would  result  from 
a  16  000  °C/s  scan.  For  this  situation,  the  diffusion  is  lim¬ 
ited  to  a  region  of  less  than  2  nm  at  the  interface.  The  dashed 
line  in  Fig.  7  shows  a  profile  calculated  for  a  120  °C/s  scan, 
where  the  predicted  mass  transport  is  similar  to  that  observed 
in  Fig.  5(b).  Significant  transport  of  Ni  has  occurred  through¬ 
out  the  20  nm  extent  depicted  in  the  model,  yet  homogeni¬ 
zation  has  not  occurred.  If  the  D(gb)  value  for  Ni  in  Ni2Si  or 
the  generalized  value  for  Ni  diffusivity  in  a  metal  matrix  are 
used,  essentially  complete  homogenization  results  for  both 


TABLE  I.  Relevant  diffusivities  at  560  °C. 


D(gb) 

Species/Matrix 

(m^/s) 

Reference 

Ni/NiSb 

3.9X10-'* 

18 

Ni/NbSi 

~5  X  10“'^ 

19 

Ni/metal  (avg.) 

5X10-''' 

20 

scan  rates.  Thus,  different  diffusivities  are  required  to  ex¬ 
plain  the  results  of  slow  versus  rapid  scan  rates.  This  sug¬ 
gests  an  active  role  of  NiSi2  at  the  slow  scan  rates  and  in¬ 
creased  influence  of  other  phases  such  as  metastable  liquid  or 


FIG.  7.  Diffusion  profile  calculated  with  D(gb)  for  Ni/NiS^.  Solid  line 
represents  calculation  for  a  16  000  °C/s  scan,  and  dashed  line  is  for  a 
120  °C/s  scan.  C/Cq  is  the  ratio  of  the  Ni  concentration  after  a  temperature 
scan  to  the  concentration  before  the  scan. 
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glass,  or  Ni-rich  silicides  (e.g.,  Ni3Si)  with  higher  diffusivi- 
ties  at  the  rapid  scan  rates.  We  know  from  SARD  and  XRD 
that  NiSi2  and  Ni3Si  are  major  product  phases  for  both  scan 
rates.  The  compositional  layering  evident  in  Figs.  5(d)  and 
5(e)  could  be  explained  by  the  influence  of  an  NiSi2-rich 
layer,  which  would  serve  as  an  effective  diffusion  barrier, 
preventing  further  homogenization  in  the  slow  scans.  The 
lack  of  layering  in  Fig.  5(a),  together  with  the  diffraction 
data  implies  interspersed  Ni3Si  and  NiSi2  in  the  rapid  scan 
products.  It  seems  unlikely  on  micro  structural  grounds  that  a 
completely  interspersed,  relatively  uniform  texture  would  re¬ 
sult  from  interfacial  mass  transport  in  the  solid  state,  since 
there  has  to  be  a  continuous  directional  component  to  the 
mass  flux  and  associated  reaction.  An  efficient  homogeniza¬ 
tion  process  seems  to  be  required,  for  which  melting  pro¬ 
vides  the  most  probable  explanation. 

D.  Effect  of  scan  rate  on  reaction  energetics 

It  is  useful  to  compare  the  observed  enthalpies  (Fig.  4) 
with  the  values  computed^  ^  for  complete  reaction  to  Ni2Si 
and  NiSi  [the  equilibrium  products  predicted  by  Fig.  6(a)  for 
0.35  mol  fraction  Si]  versus  reaction  to  the  observed  prod¬ 
ucts  (Ni3Si  and  NiSi2).  Our  calculations  suggest  8.9  /xJ  for 
the  equilibrium  reaction  and  8.3  /xJ  for  our  observed  prod¬ 
ucts.  The  maximum  value  observed  at  the  most  rapid  scan 
rates  in  Fig.  4,  i.e.,  9.8  /xJ,  is  higher  than  the  maximum 
predicted  for  the  equilibrium  phase  assemblage.  We  attribute 
this  to  added  energies  associated  with  the  amorphous  to  crys¬ 
talline  transition  in  Si,  and  to  spontaneous  grain  growth  in 
Ni.  With  these  caviats,  the  DSC  data  appear  to  be  in  reason¬ 
able  agreement  (within  15%)  of  predicted  values. 

The  fine  structure  of  peak  1  revealed  in  Fig.  3  suggests 
the  need  for  further  interpretation.  For  example,  it  is  possible 
that  a  small  amount  of  melting  (the  first  maximum  of  the 
doublet  of  peak  1)  occurs  even  in  the  slow  scans,  but  that  this 
melting  gives  way  to  crystallization  (e.g.,  NiSi2)  (the  second 
maximum  of  the  doublet),  that  impedes  further  melting  and 
homogenization.  As  the  scan  rate  increases,  the  time  avail¬ 
able  for  NiSi2  to  form  is  less,  so  that  the  amount  of  melting 
increases,  resulting  in  a  diminished  doublet.  At  the  highest 
scan  rate,  melting  and  complete  homogenization  occur  si¬ 
multaneously  with  the  crystallization  of  Ni3Si  and  NiSi2 
products,  resulting  in  a  combination  of  peaks  1  and  2. 

V.  CONCLUSIONS 

In  summary,  we  favor  the  application  of  a  metastable 
melting  hypothesis  to  explain  the  enhanced  mass  transport 
and  microstructure  produced  by  the  ultrarapid  heating  of  the 
Ni/Si  multilayers  (DSC  scan  rate  16  000  °C/s),  whereas  we 
postulate  that  solid  state  diffusion  is  responsible  for  the  mi¬ 
crostructure  produced  by  the  slow  scan  rate  (120  °C/s).  Our 


observations  and  conclusions  regarding  the  reactions  to  form 
the  silicide  phases  Ni3Si  and  NiSi2  as  final  reaction  products 
appear  to  be  consistent  with  the  DSC  enthalpies,  and  the 
relatively  low  grain  boundary  diffusivity  of  NiSi2.  These  re¬ 
sults  point  to  the  possibility  of  alternative  routes  for  practical 
processing  of  semiconductor  silicide  interconnects  under 
conditions  of  flash  annealing.  Further  combined 
microstructural/nano-DSC  experiments  will  provide  more 
details  on  the  potential  use  of  ultrarapid  heating  in  Si-based 
device  fabrication. 
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